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Abstract. A simple method of the vertical muon energy spectrum simulations have been suggested. These 
calculations have been carried out in terms of various models of hadronic interactions. The most energetic n ± - 
mesons and K*-mesons produced in hadron interactions contribute mainly in to this energy spectrum of muons 
due to the very steep energy spectrum of the primary particles. So, some constraints on the hadronic interaction 
models may be set from a comparison of calculated results with the cosmic data on the vertical muon energy 
spectrum. This comparison showed that the most energetic secondary particles production is too high in case 
of the QGSJET 11-04 model and rather low in case of the QGSJET 11-03 model. These conclusion have been 
supported by the LHC data. 


1 Introduction 

The longitudinal development of extensive air showers 
(EAS), in particular, the depth X max of its maximum de¬ 
termined by the rate of fragmentation of energy Eq of the 
primary particle. This rate depends on the interaction cross 
sections of shower particles, and on the energy spectra of 
secondary particles are generated in interactions. Obvi¬ 
ously, if the probability of particle production with ener¬ 
gies close to the energy of the incident particle is high 
then the development of the cascade slows down. Con¬ 
versely, in the case of rapid fragmentation of the incident 
energy cascade develop rather rapidly. The depth X max 
of shower maximum in many studies is the main param¬ 
eter for determining the composition of primary cosmic 
radiation (PCR) at ultra-high energies. It should also be 
noted that in the case of the slow rate of development cas¬ 
cade and hence large values of depth X max lateral distri¬ 
bution of the shower particles at the observation level be¬ 
comes narrower. Therefore, the values of signals in the 
surface and underground detectors located at large dis¬ 
tances from the shower core are decreased. It must be 
taken into account when determining the density of muons 
at large distances from the shower core and the composi¬ 
tion of the PCR found from the muon lateral distribution. 
Studies of the composition and characteristics of the en¬ 
ergy spectrum of the PCR are important components of 
theories of the origin of cosmic rays at ultra-high ener¬ 
gies. Interpretation of experimental data on the depth of 
the shower maximum X max and the observed fraction of 
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muons at a fixed distance from the shower core are earned 
out in terms of different hadronic interaction models. In 
the case of "soft" interaction generation of secondary par¬ 
ticles emitted at mostly small angles with respect to the 
projectile particle (the highest values of pseudorapidity if) 
is of importance. The models based on the Gribov-Regge 
theory [1,2] are commonly used. The dominant contribu¬ 
tion of the pomeron at ultra-high energies and other effects 
are taken into account differently in various models [3- 
5], Therefore, testing of models at the highest energies of 
secondary particles is very important for understanding of 
physics of hadron interactions and for the interpretation of 
the EAS data. At the LHC this testing is carried out in the 
experiments LHCf [6] and TOTEM [7], In cosmic rays, 
we suggest testing of models of hadron interactions with 
the help of atmospheric muon spectrum. In this case, due 
to the higher slope of the energy spectrum of the primary 
particles (a coefficient of the slope of the difference spec¬ 
trum is y = 2.75) generation of secondary particles (n ± - 
mesons and A^-mesons) with the highest energies is of 
importance. It was shown [8] that the model QGSJET II- 
03 [4] leads to a spectrum of vertical muons, the inten¬ 
sity of which is about a factor / = 1.5 times less than 
the data of collaborations L3 + Cosmic [9], MACRO [10] 
and LVD [11]. The result [8] was obtained as a solution 
of transport equations. The zr* and K ± mesons decay into 
/jE -mesons. These zri and K ± mesons are generated by 
the parent particles of several generations. Therefore, if 
some excess of the secondary particle production in sin¬ 
gular interaction of hadrons is determined by the coeffi¬ 
cient k , then, in the case of i generations, the factor / 
will be / « k'. Thus, the spectra of atmospheric muons 
are very useful tool to test the models of hadronic inter- 
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Figure 1. The primary spectrum of protons and he¬ 
lium nuclei (p+He). Dashed line - modified Gaisser- 
Honda [13] approximation. Data: solid line - AMS02 [14]; 
o ATIC 2 [15]; • CREAM [16]; ■ WCFTA [17]; 
A ARGO [18]; □ RUN JOB [19]; o TUNKA [21] (all par¬ 
ticles); ▲ SPHERE 2 [22] (all particles); x KASKADE [20] 
(all particles, QGSJET 11-03); + KASKADE [20] (all particles, 
SIBYLL 2.1). 



Figure 2. The energy spectrum of the primary protons. Dashed 
line - modified Gaisser-Honda [13] approximation. Various data 
- see the text. A AMS02 [14]; o ATIC 2 [15]; • CREAM [16]; 
□ RUN JOB [19], 


actions. However, other factors (cross sections of interac¬ 
tions, etc.) are affects result of comparison. In this article 
we propose to test the hadron interaction models by the 
very simple original method [12] with the help of observed 
atmospheric muon spectrum. 

2 Method of calculation 

The very simple original method of simulations of the en¬ 
ergy spectrum of vertical muons can be described as fol¬ 
low [12]. Let the (dI p /dE) and (dI He /dE) be the differ¬ 
ential energy spectra of the primary protons and helium 
nuclei. As for the spectrum of muons the energy per nu¬ 
cleon is of importance. So, the heavier nuclei give to this 


Figure 3. The energy spectrum of the primary helium nu¬ 
clei. Dashed line - modified Gaisser-Honda [13] approxima¬ 
tion. Various data - see the text. A AMS02 [14]; o ATIC 2 [15]; 
• CREAM [16]; □ RUN JOB [19]. 


spectra negligible contribution. In the energy range of 10 2 

- 3 • 10 6 GeV, we used the approximations ( dI p /dE) CH 
and ( dI He /dE) CH by Gaisser-Honda [13]. At energies 
above E\ = 3 • 10 6 GeV, these spectra were multiplied 
by a factor (E \ /E ) 0 5 and are refered as modified GH ap¬ 
proximation. Figure 1 shows a comparison of the sum 
of modified GH approximations [13] for the primary pro¬ 
tons and helium nuclei (dotted line) with the experimental 
data AMS02 [14] - solid line, ATIC2 [15] - hollow cir¬ 
cles (o), CREAM [16] - dark circles (•), WCFTA [17] 

- dark squares (■), ARGO [18] - hollow triangles (A), 
RUN JOB [19] - hollow squares (□). The experimental 
data KASKADE [20], interpreted in terms of the model of 
QGSJET 11-03 shows the oblique crosses (x), and in terms 
of the model SIBYLL 2.1 - straight crosses (+). The ex¬ 
perimental data TUNKA [21] for all primary particles are 
shown by hollow diamonds (o), and data SPHERE 2 [22] - 
dark triangles (A). The spectra of the primary protons and 
primary helium nuclei are presented separately in Figure 
2 and Figure 3 with various data mentioned above for a 
comparison. The notation: (A) is applied in Figure 2 and 
Figure 3 only for the data of AMS02 [14]. From com¬ 
parison with data it can be concluded that the accepted 
approximation [13] does not overestimate the flux of the 
primary protons and helium nuclei. This is important for 
the conclusions on possible uncertainties of tested mod¬ 
els. The energy spectra of vertical muons D p (E p )dE p and 
E>He(E p )dE p induced by the primary protons and helium 
nuclei are expressed by simple integrals over the energy E 
of the primary particles as follows: 

D p (E p ) • dE p = J dE-^j-S p (E p ,E)-dE p 

D He (E p ^ • dE p = JdE- 
The sum of these spectra: 

D(E p ) = (D p (E p ) + D He (E p )) 


dine 

~dE 


S He (e p , e) ■ dE p 
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Figure 4. The energy spectra of muons in showers induced by 
the primary protons with various fixed energies E (simulations 
in terms of QGSJET 11-04 model): 1 - 3,162 - 10 2 ; 2 - 10 3 ; 3 - 
10 4 ; 4 - 10 5 ; 5 - 10 6 ; 6 - 10 7 GeV. 


will be used for comparison with data [9-11]. Functions 
S P (E M ,E) ■ dEp and S He(E p , E) ■ dE fl are the differential 
energy spectrum of muons in showers induced by the pri¬ 
mary protons and helium nuclei with the fixed energy E. 
These spectra were calculated for 24 and 19 values of the 
energy E of the primary protons and helium nuclei, re¬ 
spectively, in the range of 10 2 4- 10 7 GeV. Calculations 
have been carried out in terms of two models of hadron 
interactions (QGSJET 11-03 [4] and QGSJET 11-04 [5]) 
using the package CORSIKA 7.4 [23]. The calculations 
were performed with the statistics of 10 6 events in the low 
energy region and up to 10 2 events at the highest ener¬ 
gies of the primary particles. Figure 4 and Figure 5 
show examples of the muon energy spectra S p [e p , Ej and 

SHe {Ep,E^j in the energy range of 10 2 4- 10 5 GeV calcu¬ 
lated in the terms of model QGSJET 11-04 [5] for the 6 
values of proton energies and the 5 values of helium nu¬ 
clei energies, respectively. It is seen that in the energy 
range of 10 4 4 - 10 6 GeV statistics is small. Therefore, we 
will use the energy range of 10 2 4 - 10 4 GeV for a compari¬ 
son. Dependence of the spectra S,, (Ep, Ej and S He (e p , E'j 
on the models of hadron interactions (QGSJET 11-03 (•) 
and QGSJET 11-04 (o)) are shown in Figure 6 and Figure 
7 for the primary protons and helium nuclei, respectively 
for the energy E = 10 s GeV. It is evident that the model 
QGSJET 11-04 [5] predicts the greatest density of muons, 
while the model QGSJET 11-03 [4] - the lowest one. 

3 The results of simulations 

The spectra of vertical muons D{E p ) in the energy range 
of 10 2 - 10 4 GeV for the hadron interaction models 
QGSJET 11-03 [4] and QGSJET 11-04 [5] are presented in 
Figure 8. It can be seen that the model QGSJET 11-04 [5] 
predicts the intensity of muon flux by a factor 2 4-3 higher 



Figure 5. The energy spectra of muons in showers induced by the 
primary helium nuclei with various fixed energies E (simulations 
in terms of QGSJET 11-04 model): 1 - 10 3 ; 2 - 10 4 ; 3 - 10 5 ; 4 - 
10 6 ; 5 - 10 7 GeV. 


than the intensity of the muon flux calculated in terms of 
QGSJET 11-03 [4] model. This finding is consistent with 
the results shown in Figure 6 and Figure 7. Figure 8 is also 
very clearly demonstrates an increase in steepness of the 
muon spectrum at energies E above 100 GeV. The constant 
B„ ^ 100 GeV is the decay constant for 7r-mesons in the at¬ 
mosphere. Therefore, the muon spectrum become steeper 
because 7r-mesons at energies E > B n rather interact with 
nuclei in the atmosphere than decay into muons. Compar¬ 
ison of the calculated spectra with experimental data al¬ 
lows us to test models. The ratios MC/DATA of the results 
of calculations for models [4] and [5] to the smooth ap¬ 
proximation of the data of collaborations L3 + Cosmic [9], 
MACRO [10] and LVD [11] are shown in Figure 9. It can 
be seen that these ratios are increasing from ~ 1.4 to ~ 
1.7 for the model QGSJET 11-04 [5] and are decreasing 
from ~ 0.8 to ~ 0.6 for the model QGSJET 11-03 [4] at 
muon energies in the range of 10 2 4 - 10 4 GeV. The most 
important fact is that these increase becomes higher at en¬ 
ergies Ep above 10 3 GeV for the QGSJET 11-04 model. 
For the QGSJET 11-03 model the ratios MC/DATA become 
constant at the level of ~ 0.6 at energies E fJ above 10 3 
GeV. No slowing of this increase is seen at higher ener¬ 
gies of muons. Thus, Figure 9 demonstrates a very serious 
difference between the calculated spectra and the data re¬ 
ported in [9], [10] and [11]. This difference is associated 
with a different rate of energy fragmentation of projectile 
particles in events of its interactions with nuclei in the at¬ 
mosphere. Thus, the model QGSJET 11-03 underestimate 
the probability of secondary particles production at high¬ 
est energies. The model QGSJET 11-04 overestimate this 
probability of secondary particles production at the highest 
energies. According to calculations, the main contribution 
to integrals D p yE^j ■ dE p and D He yE p J ■ dE p comes from 
secondary particles with energies in the ranges of (0.01 
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Ej, GeV 


Figure 8. The calculated energy spectra of near - vertical muons 
for various models: • QGSJET 11-03, o QGSJET 11-04. 


Figure 6. The energy spectra of muons in showers induced by the 
primary protons with the fixed energy E = 10 5 GeV calculated 
in terms of two models: • QGSJET 11-03, o QGSJET 11-04. 



Figure 7. The energy spectra of muons in showers induced by the 
primary helium nuclei with the fixed energy E = 10 5 GeV calcu¬ 
lated in terms of two models: • QGSJET 11-03, o QGSJET 11-04. 


- 0.6 )E and (0.001 - 0.1)is, where E is the energy of a 
projectile particle for the primary protons and helium nu¬ 
clei, respectively. The results of this comparison are also 
confirmed by the data of the LHCf [ 6 ] and TOTEM [7]. 
For example, the QGSJET 11-04 model [5] overestimates 
the density of charged particles dN c h/drj per unit of pseu¬ 
dorapidity at the pseudorapidity 77 = 6.345 by a factor of 
k ~ 1.3 as compared to the TOTEM data [7]. This differ¬ 
ence increases at large values 77 because of the difference 
between the slopes of the calculated curve and the data 
from [7]. The QGSJET 11-04 [5] model predicts the den¬ 
sity dNchldrj which is (18-30)% higher than the data [7] 
in the range 5.3 < 77 < 6.4. It is also important to note 
that for difference of rapidity A y - 0 , the value of the av¬ 



Figure 9. Comparison of calculated muon energy spectrum with 
data [9-11]. The ratios MC/DATA are shown: • QGSJET 11-03, 
o QGSJET 11-04. 


erage transverse momentum (pr) for 7r°-mesons is about 
50 MeV/c less than data [24]. Under the assumption that 
a similar decreasing of the (p T ) dependence is also valid 
for charged n mesons. The calculated density of muons 
at large distances from the axis shower will be underesti¬ 
mated. This underestimation was observed by the Pierre 
Auger Collaboration [25] and in Yakutsk [26]. The en¬ 
ergy spectra of photons in p-p collisions at the energy of 
\fs = 7 TeV in the pseudorapidity range 8.81 < 77 < 8.99 
are 2-r4 times above the one’s predicted by the QGSJET II- 
03 model [4]. So, all these models should be significantly 
corrected at the highest energies of secondary particles. 

4 Conclusion 

The atmospheric muon energy spectra calculated in terms 
of the QGSJET 11-04 [5] model is by a factor / = 1,7 
higher then data [9], [10], [11] at energy E fJ = 10 4 
GeV. The atmospheric muon energy spectrum calculated 
in terms of the QGSJET 11-03 [4] is by a factor / = 1,5 
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lower then data [9], [10], [11] at energy E^ = 10 4 GeV. 
So, we can conclude, these models of hadronic interac¬ 
tions should be updated at very high energies of secondary 
particles. 
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